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Abstract. Under extreme conditions of temperature and/or density, quarks and gluons are expected to 
undergo a deconfinement phase transition. While this is an ephemeral phenomenon at the ultra-relativistic 
heavy-ion collider (BNL-RHIC), quark matter may exist naturally in the dense interior of neutron stars. 
Herein, we present an appraisal of the possible phase structure of dense quark matter inside neutron stars, 
and the likelihood of its existence given the current status of neutron star observations. We conclude 
that quark matter inside neutron stars cannot be dismissed as a possibility, although recent observational 
evidence rules out most soft equations of state. 

PACS. 97.60. Jd Neutron stars - 26.60.+C Nuclear matter aspects of Neutron stars 



1 Introduction 

BNL-RHIC is engaged in a voyage of exploration and dis- 
covery in the high-temperature, low baryon density regime 

\ of QCD's phase diagram [HE]- One of the central aims of 
this program is to characterize the deconfinement and chi- 
ral phase transition in QCD at temperatures reminiscent 
of the hot Big Bang. It is widely believed that a strongly 
interacting medium of quarks and gluons, displaying ideal 
liquid-like behaviour, has been created in the most ener- 
getic central Au-Au collisions at RHIC[21I21[S]. However, 
the lifetime of this phase is of the order of lOfm/c, re- 
quiring penetrating experimental probes that carry the 
imprint of the early hot partonic phase, and which are 

' not washed out by the hadronization process [BUTUH]. 

The dense interior of neutron stars provides a complemen- 
tary testing ground for quark deconfinement. The central 
densities inside neutron stars can be as high as 5-lOpo 
(po = 2.5 x 10 14 g/cc is the nuclear saturation density), 
and nucleons overlap to an extent that quarks and glu- 
ons become the effective degrees of freedom. Under such 
extreme conditions of density, it is possible that strange 
quark matter is energetically more stable than nuclear 
matter [9KT0] . If this is the case, there is a critical pres- 
sure at which a first order phase transition from nuclear 
to quark matter will occur. Quark matter can then com- 
prise an arbitrary fraction of the star, from zero for a pure 
neutron star to one for a pure quark star, depending on 
the equation of state of matter at high density. 

In these proceedings, we outline the rationale behind the 
possible stability of three flavor (up, down and strange) 
quark matter at high baryon density and review results 
from strangelet searches in heavy-ion collisions from AGS 
to CERN-SPS and RHIC. We consider the effects of finite 
size and interface energy corrections to the phase structure 



of stable strange matter and focus on the consequence 
for the surface structure of quark stars. We emphasize 
recent developments in neutron star observations that can 
shed light on the possible existence of quark matter inside 
neutron stars. 

2 Strange quark matter in heavy-ion collisions 

The rationale behind stable strange quark matter is the 
Witten hypothesis 10 , which argues that the introduction 
of strangeness in up and down quark matter reduces Pauli 
repulsion by increasing the flavor degeneracy from up and 
down to up, down and strange quarks, ensuring in the pro- 
cess also a lower charge-to-baryon ratio for strange quark 
matter compared to nuclear matter. While this hypoth- 
esis is clearly not borne out for small baryon numbers, 
where strange baryons are definitely heavier than their 
non-strange counterparts, there is no observational evi- 
dence to suggest that this is also the case in bulk strange 
quark matter. Nuclei would not spontaneously decay to 
strange matter even if the latter was more stable, since 
that would require ~ A weak reactions to occur simulta- 
neously in a nuclear volume containing A nucleons. This 
conversion would happen much more easily in the interior 
of neutron stars, where pressures and densities are supra- 
nuclear. The critical question regarding quark matter in 
neutron stars is then whether the central density of neu- 
tron stars is large enough so that strange quark matter 
becomes the ground state of strongly interacting matter. 
This question evades a precise answer because QCD is 
still not sufficiently well understood at neutron star densi- 
ties. Lattice methods fail at such high densities due to the 
complexion of the measure involved in importance sam- 
pling. In the absence of concrete results from lattice stud- 
ies of QCD at finite density and zero temperature, simple 
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model-dependent studies [TT] admit a parameter window 
(the parameters being the strange quark mass, the strong 
coupling constant and a phcnomcnological Bag constant) 
within which bulk strange quark matter is stable, even at 
zero pressure. If true, this implies that, if central densi- 
ties inside neutron stars are large enough to create two- 
flavor (up and down) quark matter, or if a small nugget of 
cosmological/cosmic-ray origin ("strangelet") enters the 
star, the entire neutron matter inside the star will con- 
vert to strange quark matter by absorbing neutrons and 
equilibrating strangeness. 

Strangelet searches in terrestrial materials, cosmic rays or 
as by-products of neutron star-neutron star collisions have 
thus far yielded negative resultsp^, 13 . Even if strange 
quark matter is stable in bulk, it may be destabilized 
by prohibitive surface and Coulomb energy costs so that 
strangelets do not survive until the present day although 
they may have existed in the hot and dense epoch of 
the early universe. If so, conditions in the forward rapid- 
ity regime of ultra-relativistic nucleus- nucleus collisions 
(and mid-rapidity at fixed target experiments) may be 
able to create strangelets for a short while before they 
evaporate [T4lll5| . The experimental signal searched for is 
a particle with a large mass-to-charge ratio that, owing 
to its large rigidity, would not be deflected by magnetic 
fields, and would be able to reach the zero-degree calorime- 
ter (ZDC). There, it would produce a shower originating 
from a single point, unlike spectator neutrons which are 
dispersed in the transverse plane due to Fermi motion. 
While the strangelet search at NA52[14] at the CERN- 
SPS was sensitive to long lived strangelets(r ~ fj,s), the 
corresponding experiments at AGS [16] and at RHIC[l5] 
were sensitive down to lifetimes of t ^50ns. From these 
experiments, the production rate of strangelets was lim- 
ited to less than one in 10 7 — 10 9 central collisions for 
strangelets exceeding a mass of 30 GeV/c 2 and lifetime 
greater than a few nanoseconds. 

Such a low probability for producing strangelets in heavy- 
ion collisions is expected on theoretical grounds as well. 
Various models have been examined as a mechanism to 
produce stable strangelets (see Ref [17] for a review) . The 
coalescence mechanism involves forming a clump of strange 
matter by overlap of a sufficient number of baryons (of 
appropriate strangeness) with small relative momentum. 
This is highly improbable at collider energies. The thermal 
model, which has its parameters, temperature and bary- 
ochemical potential, tuned to reproduce observed parti- 
cle ratios at chemical freeze-out, predicts less than one 
strangelet for every 10 27 collisions for a strangelet with 
Z/A ~ 0.1 and mass 20GeV/c 2 . This production mech- 
anism falls off rapidly with increasing collider energies. 
The QGP distillation method relies on the enhancement 
of strange quarks in a QGP followed by evaporation of 
the baryon rich QGP through nucleons, thereby distilling 
strangeness and enabling a cool and stable strangelet to 
emerge. However, baryon rich and strangeness-rich regimes 
are well separated in rapidity, and there is clear evidence 
that the QGP cools through rapid, approximately adi- 



abatic expansion rather than by evaporation|18] . Alto- 
gether, it is very unlikely that stable strangelets would 
be produced in a heavy-ion collision. 



3 Strange quark matter in neutron stars 

If strange quark matter is stable only at very high baryon 
number (A ^> 10 7 ), a neutron star with 10 57 baryons is 
a natural candidate where such matter can exist. Two 
possibilities then arise: (i) quark matter is stable in bulk 
at some large value of the pressure. In this case, a first- 
order transition is likely to occur at some depth (density) 
inside the neutron star and quark matter is admixed with 
hadronic matter in a mixed phase whose structural details 
are determined by surface and Coulomb effects, (ii) quark 
matter is stable in bulk even at zero pressure (still at finite 
density, since it is self-bound). This would imply that all 
neutron stars are really strange quark stars, with possibly 
a thin layer of hadronic matter at surface. Let us examine 
these possibilities in more detail below: 

(i) In the event that a first order phase transition oc- 
curs inside the star, a mixed phase of nuclear and quark 
matter can occupy a significant portion of the star's inte- 
rior. This conclusion follows from the fact that there are 
two conserved quantities, electric charge and baryon num- 
ber, which can be arranged differently in the two phases, 
quark and nuclear, at different equilibrium pressures. Thus, 
we expect a gradually increasing proportion of quark mat- 
ter with increasing depth inside the neutron star. The 
structure and size of the rarer phase (droplets/rods/slabs) 
at a given density depends on the surface tension between 
the two phases, the curvature energy and the smallest De- 
bye screening length. While positive surface tension and 
curvature energy tend to disfavor small sizes, Coulomb 
energy disfavors large sizes, leading to deformed struc- 
tures when Debye screening effects are included. If sur- 
face tension and Coulomb costs are prohibitively large in 
the mixed phase, the standard picture of a sharp interface 
with a density discontinuity between hadronic and quark 
matter induced by gravity, is applicable, even though it is 
not the minimum energy configuration in bulk matter. 

(ii) If strange quark matter is stable in bulk even at 
zero pressure, what we call neutron stars are really quark 
stars that contain quark matter almost upto the surface. 
At scales where the strange quark mass to 2 /4/iq <C 1, 
with /iQ the quark chemical potential, quark matter is 
effectively neutral with equal numbers of up, down and 
strange quarks. Near the surface of the star, however, 
where m 2 / is not small, electrons are required to make 
up the deficit of strange quarks in order to form a neu- 
tral object. Microscopically, the electron distribution at 
the surface is governed by electrostatics on the length 

scale l e ~ 1/\J apl^ 3 ~ lOOOfm (p e is the number den- 
sity of electrons and a is the fine structure constant), 
while quarks are bound by the strong force (QCD scale 
~ lfm). Consequently, charge neutrality at the surface 
is impossible at scales smaller than a lOOOfm in a pic- 
ture where quark matter is assumed to be homogeneous. 
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The electrons in this case distribute themselves accord- 
ing to the laws of electrostatics and mechanical equilib- 
rium. They form a thin charged skin atop the star, which 
is held to the surface by enormous electric fields E ~ 
10 16 V/cm[]15]. Such an enormous electric field is expected 
to emit electron-positron pairs at sufficiently high temper- 
atures T > 10 10 K, thereby producing a dramatic signal of 
hot quark stars. This signal is also transient, with a life- 
time of a few days, since the star cools rapidly to lower 
temperatures, shutting off the pair emission |20j. Thus, ob- 
servation of this signal is highly improbable from an as- 
trophysical viewpoint. 

There is a more appealing alternative for the surface struc- 
ture which takes into account the fact that a two-phase 
system can be globally rather than locally charge neutral [21] 
Relaxing the condition of local charge neutrality allows to 
reduce the strangeness fraction in quark matter at small 
HQ-, thereby lowering its free energy. Global charge neu- 
trality is achieved in a mixed phase by having the phase 
fractions vary as a function of the pressure. This hetero- 
geneous phase is favored when surface and Coulomb en- 
ergies are negligible, as shown in the model independent 
approach followed in ref[22]- Since He <§; HQ for all rea- 
sonable equations of state describing dense quark matter, 
the quark pressure may be expanded in powers of He/l^Q- 
To second order in /i e //iQ, it is given by 



P(pQ, Ve) = Po(fJ-Q) ~ n Q {fI Q ) He + ^ Xq(^q) £ , W 

where tiq{hq) — —dP/dfi e is the positive charge den- 
sity, Xq{^q) = d 2 P/dnt is the charge susceptibility and 
Po is the pressure of the electron-free quark phase. They 
depend on fiQ, "*s, and strong interactions. Typical val- 



ues in the Bag model description are hq 



Xq = and P a = 3(hq - m 2 s H Q )/^ - B, where 

B is the bag constant. At fixed \iQi from Eq. [1] the quark 
pressure P q is zero and quark matter is positively charged 
when /i e takes on the value 



XQ 



(1 — yl — £) where £ 



Q 



A mixed phase is possible when < £ < 1. In this regime, 
the mixed phase has lower free energy (larger pressure) 
than homogeneous matter. The mixed phase is however 
penalized by Coulomb, surface, and other finite size con- 
tributions to the energy. At zero temperature and pres- 
sure, the magnitude of the change in Gibbs free energy 
per quark in going from a homogeneous to a mixed phase 
should be more than the surface tension, i.e. the energy 
cost of creating a droplet surface, Then, the mixed phase 
is preferred over homogeneous quark matter. This critical 
surface tension is [23] 



C crit 



0.8n Q 



12 



3/2 
TTa X Q 



(3) 




Fig. 1. Surface plus Coulomb energy cost as a function of 
nugget size. The optimum size of the quark nugget for the 
choice (ii) of parameters described in the text is 8.25fm.The 
inset shows the quark and electron number densities, as well 
as the electrostatic potential cj> inside the nugget. 



In the context of the Bag model for dense quark matter, 
the condition for forming a mixed phase becomes 



a < 12 



150 MeV 



^ MeV/fm 2 . 



(4) 



Using estimates of the surface energy of strangeletspMI 
[25]: (i) a ~ 8 MeV/fm 2 for m s = 150 MeV and hq ~ 300 
MeV; and (ii) a ~ 5 MeV/fm 2 for m s = 200 MeV at 
Hq ~ 300 MeV. The condition in Eq. [4] implies that a ho- 
mogeneous phase is marginally favored for m s = 150 MeV 
while the structured mixed phase is favored for m s = 200 
MeV. The sensitivity to m s in Eq. [4] and uncertainty in 
other finite size effects can alter these quantitative esti- 
mates. If the structured phase is favored, it will be com- 
posed of quark nuggets immersed in a sea of electrons. The 
size of the quark nuggets in this phase is determined by 
minimizing the surface, Coulomb and other finite size con- 
tributions to the energy (see FigJT]). At low temperature, 
this mixed phase will be a solid with electrons contribut- 
ing to the pressure while quarks contribute to the energy 
density - much like the mixed phase with electrons and 
nuclei in the crust of a conventional neutron star. This 
modified picture of the strange star surface has a much 
reduced density gradient at surface and negligible electric 
field unlike the old paradigm for quark stars. In the mod- 
ern viewpoint, there is no need for the electron skin or 
associated large electric fields, since matter at the surface 
is globally neutral. The observed luminosity from such a 
surface will be very different than from a charged surface 
with an electron skin [26]. 

It is also of interest to estimate the radial extent AR of 
the mixed phase crust. This is because some neutron stars 
exhibit "glitches" in their rotation, when they suddenly 
spin-up before gradually resuming spin down. If strange 
stars have a large crust, where a superfluid and a lattice 
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structure can co-exist, they could explain this glitching 
mechanism. Using Newtonian approximations to Hydro- 
statics, the size of the crust is given bv[2"2] 

AR = f-^-R, (5) 
Rs XQ e o 

where R s = 2GM /c 2 ~ 3(M/M S ) km is the Schwarzschild 
radius of the star, R is the radius of the star and eo is 
the energy density inside a quark nugget. For m s = 150 
MeV and fi Q ~ 300 MeV, the Bag model with B = 65 
MeV/fm 3 yields n Q ~ 0.045 mi" 3 , \Q - 46 MeV/fm, and 
e ~ 283 MeV/fm 3 . From Eq.©, Z\i? ~ 100 meters for a 
star with mass M = 1.4M . The Newtonian estimate for 
AR is close to a more accurate value for AR obtained by 
solving general relativistic equations for hydrostatic equi- 
librium numerically [52"]. Fig. [5] shows the density profile of 
the crust thus obtained. 

The interior structure of a quark star can be quite com- 
plicated. Recently, a lot of progress has been made in un- 
derstanding QCD at asymptotically high densities. In that 
regime, where pertubative studies are reliable, quark mat- 
ter is believed to be in a color-flavor-locked (CFL) phase, 
characterized by quark pairing with a completely gapped 
spectrum. Such a phase is an electromagnetic insulator in 
bulk and admits no electrons, even when stressed by small 
quark masses [27.. If dense quark matter indeed exists in- 
side neutron stars, where densities are well above nuclear 
matter density but below the density where perturbative 
QCD is expected to be valid, the ground state of quark 
matter is uncertain. Nevertheless, in such a "hybrid star" , 
attractive interactions between quarks will lead to the for- 
mation of a color superconducting state, characterized by 
quark pairing and superfluidity. The singlet pairing gaps 
can be as large as 100 MeV [55] and they modify transport 
properties due to the presence of collective excitations be- 
low the scale of the gap. At densities relevant to neutron 
stars, with hq ~ 500 MeV or less, and with the physi- 
cal requirements of charge and color neutrality, the pair- 



ing pattern can be quite complex involving phases with 
gapless modes for certain quark quasiparticles. Of partic- 
ular interest are the crystalline phases [29], where quarks 
with different Fermi surfaces pair at non-zero momentum, 
resulting in an inhomogenous but spatially periodic or- 
der parameter. The crystalline structure may also serve 
as sites for pinning rotational vortices formed in the su- 
perfluid as a result of stellar rotation, and could generate 
the observed glitch phenomena in neutron star spin-down. 
The gapless and crystalline phases can also lead to temper- 
ature dependences that are modified from the usual forms 
in ungapped quark matter. These phases also have unique 
dispersion relations for certain quark quasiparticles, and 
consequently, a specific heat per unit volume that is also 
different from ungapped quark matter. These two factors 
imply a change in the stellar cooling curve that can be 
confronted by observations [30 . The question as to which 
is the preferred stable state of quark matter at intermedi- 
ate densities and physical strange quark mass remains an 
open one at this time. 



4 Constraints from neutron star observations 

Neutron star observations can help in constraining the 
equation of state of dense matter, and also in distinguish- 
ing between different models for the crust as discussed 
above. Individual neutron star masses are most precisely 
determined by measuring post-Keplerian orbital param- 
eters in close binary systems. Neutron star masses thus 
determined lie in the range 1.18-1.44 M and have errors 
of less than a tenth of a percent. If the radius can also be 
ascertained to high accuracy, the equation of state of dense 
neutron star matter can be pinned down. In practice, there 
are several complications that make radius measurements 
a challenge so that it is only possible to infer the radius 
at infinity which is related to the true radius of the star 
through the relation R^ = (1 + z) R with the red-shift 
factor (l + z) -1 = yl — 2GM/Rc 2 . Consequently, instead 
of measuring a radius we can only infer a relation between 
mass and radius. It has recently been realized that com- 
pact objects in low-mass X-ray binaries (LMXBs) which 
exhibit X-ray bursting behavior may provide a promising 
new avenue to determine, simultaneously, both the mass 
and radius of a neutron star [31], thereby setting firm con- 
straints on the equation of state of dense matter. In these 
objects, there is the potential to observe, in addition to 
the quiescent luminosity that can be used to infer Roo, 
the Eddington luminosity during the burst and the grav- 
itational red-shift through direct observation of the shift 
in identifiable atomic absorption lines in the atmosphere. 
The simultaneous determination of mass (2.10± 0.28 M©) 
and radius (13.8±1.8 km) [31] of the bursting LMXB EXO 
0748-676 eliminates most soft equations of state and is 
compatible only with the stiffest neutron or quark matter 
equations of state. It does not rule out all quark matter 
equations of state, therefore, hybrid stars or strange stars 
remain viable [55], However, this and other heavy neutron 
star candidates and the rather large inferred radius for 
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EXO 0748-676 disfavors the scenario in which significant 
softening due to a phase transition at high density occurs. 
These recent developments are reviewed in [33ll34| . 



5 Conclusions 

The existence of stable strange quark matter remains a dif- 
ficult proposition to veto or verify. The true ground state 
of strongly interacting matter at high density is as yet 
unknown, but useful constraints are placed by recent neu- 
tron star observations. In addition, observed astrophysical 
phenomena such as glitches, quasi-periodic oscillations in 
accreting neutron stars, thermal radiation from quiescent 
LMXBs and seismic vibrations during magnetar flares can 
potentially yield valuable information about the neutron 
star crust and interior, thereby constraining the equation 
of state at high density. It is apparent that recent neutron 
star observations disfavor the appearance of soft exotic 
(non-hadronic) matter at high density. They do not as yet 
completely rule out all plausible quark matter equations 
of state. Given that strangelet searches in astrophysical 
as well as terrestrial (heavy-ion collider) experiments have 
yielded null results, it appears that our access to the un- 
derlying degrees of freedom in QCD may be limited to 
the transient partonic fireball created in ultra-relativistic 
heavy-ion collisions. 
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Andrew Steiner, Rachid Ouyed, Andreas Schmitt, Mark 
Alford and Giorgio Torrieri; and I thank the organizers of 
Hot Quarks 2006 for an invigorating workshop. This work 
is supported by US DOE grant DE-FG02-93ER40756. 
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